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METHODS AND DEVICES FOR GROUPING CELLS 

BACKGROUND OF THE INVENTION 
[0001] As the name implies, cellular networks are made up of many separate cells. 
Commonly, cells are grouped into mobile switching center (MSC) domains in order to 
route calls between individual cells and the remainder of a wireless network (e.g., 
cellular, Personal Communications System (PCS), etc.). When a user of a wireless 
5 device (e.g., cell phone, pager or PDA) moves from one location to another, the network 
must track her location, or more precisely, the location of her mobile device. Typically, 
there are thousands of such movements each second in a modern wireless network. A 
wireless network tracks the number of such movements and associates a network cost to 
them. This cost is referred to as an updating cost. 

10 [0002] While updating costs are associated with tracking the movement of users, 
other costs, referred to as paging costs, are associated with finding users. For example, 
when one user sends a call from one wireless network to a user within another wireless 
network, the receiving network must identify and locate the intended recipient of the call 
within its many cells. Typically, existing wireless networks send a page to each cell 

15 within an MSC where the recipient might be located. For example, if there are 100 cells 
within an MSC, a typical wireless network will send a page to all 100 cells even though 
the recipient is located in only one of the cells. The network tracks these pages and 
assigns a paging cost to each page. 

[0003] Owners, operators and others involved in wireless networks desire to 
20 minimize or reduce the updating and paging costs associated with their networks. To do 
so, it is desirable to devise techniques which partition or group cells in such a way that 
both updating and paging costs are minimized. This has been a daunting challenge for it 
has been thought that even the mathematical statements (i.e., equations), which represent 
the problem have, heretofore, only been hypothetically solvable using exhaustive 
25 searches/iterations (none have actually been completed). 

[0004] Finding a way to partition cells is a "balancing act" of sorts. Ideally, updating 
costs could be minimized by limiting the number of so-called "location area boundaries" 
that a user traverses when she moves from place to place. One way of doing this is to 
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include every cell in one large "location area" (LA). As one might surmise, while this 
minimizes updating costs, it greatly increases paging costs because the larger the number 
of cells, the larger the number of pages that must be sent. The reverse is also true (i.e., 
creating new location areas by removing cells from an original location area may 
5 decrease paging costs but increase updating costs). 

[0005] To date, though some have been able to formulate mathematical equations 
representing the parameters involved in partitioning cells, none have been able to 
guarantee that their solutions are efficient. 

SUMMARY OF THE INVENTION 

10 [0006] The present invention uses a linear program which, after its solutions are 
rounded, provides an efficient means of partitioning or grouping cells in order to 
approximate an optimum sum of both updating and paging costs for a wireless network 
or the like. 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 [0007] FIG. 1 depicts a simplified view of cells within a wireless network. 

[0008] FIG. 2 depicts a simplified view of cells grouped into location areas. 

DETAILED DESCRIPTION OF THE INVENTION 

[0009] Referring to FIG. 1, there is shown a wireless network 1 which comprises 
cells a-j (where "j" is the last cell in the network). The task at hand is to partition the 
20 cells a-j such that the updating and paging costs associated with the operation of wireless 
network 1 are minimized. 

[0010] FIG. 1 also shows a second wireless network 200 comprising cell aa. Greatly 
simplified, the present invention seeks to minimize the sum of the paging costs (which 
occur, for example, when a user within cell aa attempts to contact a user within wireless 
25 network 1) and updating costs (which occur, for example, when a user within network 1 
moves from one cell (e.g., cell a) to another cell (e.g., cell b) within network 1). 
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[0011] In one embodiment of the present invention, cells a-j are grouped, clustered 
or partitioned by: generating a linear program representing a sum of a contribution of 
each cell a-j to a total paging cost and a total updating cost taking into account some 
constraints; and assigning each cell to a group (e.g., location area) in accordance with 
5 solutions of the linear program. 

[0012] Having just provided an overview of the present invention, a more detailed 
explanation will now be presented. 

[0013] Wireless network 1 may be represented by a graph G(V, £), where Vis the set 
of [V]ertices or nodes of the graph and E is the set of [£]dges. Each vertex represents a 
10 cell in the network, and there is an edge or boundary between vertices in the graph if the 
cells in the network are adjacent. The number of vertices (represented by |v| ) equals n 

and the number of edges \E\ equals m. Each vertex has a unique identification number 

ie [l ... n] and a weighted value, w,-, associated with each cell that reflects each cell's user 
population (e.g., the number of users which are operating within a given cell during a 
15 given period of time, such as rush hour). Every edge, denoted (/, j) e E , has a weighted 
value, nj, associated with each edge between adjacent cells and which specifies the 
amount of user traffic between endpoints of each edge in both directions over a given 
period of time. For completeness sake, nj is set equal to 0 for every pair of nodes 
(iJ)eV«,jM E . 

20 [0014] It follows that the weighted value, w,-, is representative of the paging cost 
associated with a cell because any time a call comes in to a user in that cell, some paging 
must be carried out. Similarly, the weighted value, r,j, represents update costs if cells i 
and j are in different location areas. 

[0015] FIG. 2 depicts a greatly simplified illustration of how weighted values w, and 
25 Tij are used. Cells are shown in FIG. 2 grouped into location areas A and B. Location 
area A is given a weight w A equal to 8. Again, greatly simplified, this indicates that there 
are 8 cells within location area A. Location area B is given a weight wb equal to 4 . 
Suppose further that the number of times a user traverses the boundary or edge between 
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location area A and adjacent location area B is equal to 5 for a given time period, then 
tab - 5. 

[0016] One way of understanding the values shown in FIG. 2, is to associate each w/ 
with the number of cells within a location area which may have to be paged. Because 
5 there is a greater number of cells within location area A, intuitively it may occur to the 
reader that the paging costs associated with location area A is higher than location area B; 
that is correct. 

[0017] On the other hand, user traffic weight r i} represents the number of times a user 
passes from one adjacent location area into another. For example, if a user passes from 

10 location area A into location area B that amounts to one unit of traffic. Five units of 
traffic between location area A and location area B indicates that in average five users 
traverse the edge over a given time period. Perhaps more commonly, when r AB = 5 this 
represents five separate traversals of the same edge by five separate users. It should be 
understood that a user may traverse an edge in either direction, e.g., three users travel 

15 from A to B and two users travel from B to A, the total adding up to five units of traffic 
between the endpoints of the edge represented by the in FIG. 2. 

[0018] The first step in arriving at a solution which has some degree of reliability is 
to first formulate an expression of the problem which is solvable without the need to 
resort to exhaustive searches. To that we now turn. 

20 [0019] Generally, the paging cost can be defined as: 

Page _ Cost(L) = A • C P ■ £ \s\ • £ wi (1) 

where L is equal to {Sy, S*}, X denotes a user incoming call rate, and C p indicates the 
cost of paging a single cell. 

[0020] In Equation (1), the paging cost of a single location area having |s| cells is the 
25 product of the incoming call rate A • X/ € s wi times the cost of paging all of the cells 
within every location area, C P • \s\ . 


Atty Docket No. 29250-000998/US 
Lucent Reference: 125670/Bejerano 4-1-1-2 

Page 5 

[0021] Similarly, an updating cost can be defined by: 

Update _ Cost (L)=-CV^ £ ( 2 ) 

where C M is the cost of a single update operation and C u times nj is the update cost caused 
by traffic between cells i and j. 

5 [0022] In Equation (2), the total updating cost within a wireless network, such as 
network 1, is the amount of traffic between location areas times C u . 

[0023] To minimize the overall costs of a wireless network, the sum of the paging 
costs, represented by Equation (1) and updating costs represented by Equation (2), must 
be minimized. This can be represented by: 

10 Cost(Lopr) = min l {Update _ Cost(L) -h Page _ Cost(L) } (3) 

As will be understood by those skilled in the art, Equation (1) is a non-linear equation 
because it contains the product of two variables. The two variables being the size of the 
location area \s\ and the cell weight w t . Consequently, Equation (3) is also non-linear 

since it contains Equation (1). To formulate a solvable expression it is necessary to 
15 convert the non-linear program given by Equation (3) into a linear program. 

[0024] In one embodiment of the present invention, a linear program is generated 
from the non-linear program given by Equation (3) in order to provide approximations 
for the sum of the paging and updating costs. The details of how the present inventors 
generated a linear program from a non-linear one will be presented in more detail below. 
20 Before continuing, however, some additional comments are noteworthy. 

[0025] In actual wireless networks, geographical considerations and network 
infrastructure may impose certain size (i.e., size of a location area) and connectivity 
constraints on the network. These considerations must be taken into account in order to 
generate a meaningful formulation. For example, because all the cells of a location area 
25 are connected to a single MSC (such as MSC 2A, 2B . . . 2N shown in FIG. 1), neither the 
size nor total population of a location area should exceed the capacity of an MSC. 
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[0026] These size constraints are taken into consideration by introducing two bounds, 
K niax and W niax , on the maximal cell number and maximal population size of a location 
area, respectively. For every SeL , it is required that |s|<ATmax and 

y^i e s Wi< Wmax . In another embodiment of the present invention, more general 
5 constraints bound, for each vertex i e V , the size and weight of a location area 
containing /, i.e., \LA(iJ < Ki and £ ; e law wj < Wi . 

[0027] Topological considerations also place constraints on how cells are grouped or 
clustered into location areas. For example, it may be that certain cells must reside in the 
same location area or other cells must reside in separate location areas. These constraints 
10 are defined by constants, by, for every pair of cells i, j e V , namely: 


bij = 


1 If i and j must be in different LAs; 
- 1 If i and j must be in the same LAs; and 
0 Otherwise. 


[0028] The constants by can be represented by a matrix B = {by}. Any location area 
plan (LAP), denoted by the symbol L, that satisfies both size and connectivity constraints 
is called a feasible location area plan. 

15 [0029] In more detail, a feasible LAP can be generated as follows. Given a graph 
G(V,E), weights w/ and ry for every node is V , edge i, j e V , bounds K fnax , W niax , and 
connectivity matrix B, an L can be generated such that, 


Cost(L) = min 


2 , K } (4a) 


A'Cp'^Se i,(|S|-Zie S Wi)\ 


subject to: 
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VSuSieL: Si^Si = 0 

VieV: l<\LA(i)\<K max 

Vie V: Z;ew(o wj<W 


max 


(4b) 
(4c) 


(4d) 


5 


ViJeV,bij = l: 
ViJeV,bv = -l: 


LA{i) * LA(j) 
LA(i) = LA(j) 


(4e) 
(40 


[0030] In sum, Equations (4a)-(4f) represent a non-linear program which takes into 
consideration the constraints discussed immediately above. It is Equations (4a-4f) which 
must be converted into a linear program. Before continuing, it should be noted that the 
term "linear program" is a term of art known by those skilled in the art of partitioning 
10 cells. Though a software program may eventually be written to carry out the linear 
program described herein, the two are not synonymous. The present invention seeks to 
minimize the sum of updating and paging costs using a linear program. The techniques 
set forth herein to minimize these costs are fundamentally different from existing 
techniques. 

15 [0031] In one embodiment of the present invention, the next step in generating a 
linear program begins with defining a pair (V, d), where V is a set and d is a non- 
negative function d:V R x V — ► R (called a semi-metric) if and only if d satisfies the 
following three conditions: (i) d^ = dy x for all iJeV (symmetry); (ii) da = 0 for all 
ie V ; and (iii) d{ < d t k + d# for all i, j,k e V (triangle inequality). 

20 [0032] As provided by the present invention, the function d:V R x V — > R is used to 
partition a graph representing multiple nodes into location areas. For every pair of nodes, 
i,jeV,a variable dij e {0, l} is defined, such that, 


[0033] Those skilled in the art will recognize that conditions (/) and (ii) from above 
25 are met. Condition (iii), which requires triangle inequality, is also met as follows. 
Consider an assignment to the variables dy that induces a semi-metric. This assignment 
defines a partition of the graph into location areas in a natural way. For every node ieV , 



a first value, e.g., 1, if i and j belong to different LAs; or 
a second value, e.g., 0, if / and j belong to the same LAs. 
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the location area containing it, LA(i), is defined by, LA(i) = {j\js V a dy = 0}, i.e., all 

nodes that are zero distance from node /. Edges (/, j) e E for which dy = 1 are referred to 
as cut edges and L denotes the partitioning of LAs induced by variables dy. 

[0034] Recall the matrix B defined above. It introduces constraints to ensure that a 
5 program generates meaningful results (i.e., practical groupings of location areas). In a 
further embodiment of the present invention, for every pair /, ye V the constraints by < 
dij + 1 are applied. If nodes i, j should be in the same location area then by = -1, and the 
non-negativity constraint on dy yields dy = by +1=0. Similarly, if nodes i, j are required 
to be in different location areas then by = 1. Because dy < 1, it follows that dy = by = 1. 
10 If nodes /, j are not constrained then by = 0, and dij e {0, l} . 

[0035] In one embodiment of the present invention, the paging cost of a wireless 
network in terms of the variable dycan be defined as follows: 

Page _ Cost(L) =A-C P - ^(l-rfy)- wj (5) 

where it will be recognized by those skilled in the art that node i is paged whenever there 
15 is an incoming call to a user within the location area the node is assigned to, LA(i). 
Therefore, the network involved performs A • ^(l-cfy)- wj paging operations over a 

given time interval. 

[0036] The updating cost can also be represented in terms of the variable dy as 
follows: 

20 Update _ Cost(L) = - • C« • ^ dij - ■ m (6) 

^ ijev 

where the updating cost(s) shown in Equation (6) is simply the cost of the partitioned 
edges multiplied by the cost of a single updating operation. 
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[0037] Combining Equations (5) and (6) we arrive at a linear program substitute for 
the non-linear program given by Equation (3), which represents the total cost we are 
trying to minimize or optimize, namely, 

min<U-CV y £ j (l-dv)wj+--Cu- ^durA (7) 

5 [0038] In further embodiments of the present invention, before a final linear program 
can be formulated, size and connectivity constraints must be applied to Equation (7). 
Recall the bounds and W jnax defined above. In this instance, for each node ie V, 

LA(i) = {j\j E V a dij — 0}, \LA(i)\ = £ (1 - dij) (number of nodes in LA(/)), and 
w (LA(0) = ^ (1 - dij) • wj (weight of nodes in LA(0), we can enforce the size 
10 constraints on each location area by adding for each node ie V the constraints, 
|LA(0| < Kmax and w(LA(i)) < W max . Equation (7) can be rewritten as follows: 


min^yl-CV ^(l-dij)wj + — 'Cw ^dyrA (8a) 
subject to: 

15 ViJ,keV: dij + djk>dik (8b) 

Vi f j e V : bu < dij < bo + 1 (8c) 

VieVl ^ (1 - rfy ) < K max (8d) 

Vi g V : £ (1 - rfy) w* < W max (8e) 

ViJeV: dos {0,1} (8f) 


20 [0039] In a further embodiment of the present invention, Equations (8a)-(8f) 
represent a linear program which can eventually be used to generate (i.e., partition) cell 
clusters or groups that can be used to derive efficient approximations of the sum of 
paging and updating costs. More specifically, given the restriction placed on dy by 
Equation (8f) (i.e., that it must be a 1 or 0), Equations (8a-8f) can be said to represent a 

25 linear, integer program. 


[0040] Before continuing, it should be noted that the linear, integer program given by 
Equations (8a)-(8f) may be used in applications other than wireless communications. 
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Any application which requires the partitioning of a graph (as that term is known by 
those in the art) using at least two parameters (e.g., cell edges/boundaries and cell 
size/properties) may make use of Equations (8a)-(8f). 

[0041] Though the formulation of a linear, integer program is a novel 
5 accomplishment all by itself, more is needed before solutions can be generated and cells 
partitioned/grouped. Equations (8a)-(8f) are still considered very difficult to solve (so- 
called "NP hard")- More specifically, they would still require exhaustive searches or 
iterations to solve. In a further embodiment of the present invention, the restrictions 
placed on the linear integer program by Equation (8f) may be relaxed to allow 
10 approximate solutions to be generated in polynomial time (i.e., without exhaustive 
searches). 

[0042] In a further embodiment of the present invention, the restrictions placed on 
the linear program by Equation (8f) may be relaxed to allow for fractional solutions, i.e., 
we only require that for all V, dije [0, l] . The linear program generated by the 
15 present invention contains only O(n) variables and O(n) constraints. As will be 
recognized by one skilled in the art, the value of an fractional solution is a lower bound 
on the value of an optimal integral solution. 

[0043] Though the generation of fractional solutions allows the linear program 
represented by Equations (8a)-(8f) to be solved without the need for an exhaustive search, 
20 once the fractional values are generated they must be used to partition cells or more 
specifically, to assign a cell to one or more location areas. 

[0044] In yet a further embodiment of the present invention, the fractional values are 
rounded into integer values (e.g., 1 or 0) in order to so partition or assign each cell. 

[0045] In one embodiment of the present invention, one of many techniques which 
25 may be used to round fractional values into integers is referred to as a region growing 
(sometimes called "ball growing") technique. 

[0046] Some notation useful in understanding such a technique is as follows. A ball 
b(i, r) of radius r, i.e., a subgraph induced by these nodes, and a function (r - dij)/djk 
having edges (/, k) with only one endpoint jeb(i,r) (note that d tj is defined for all 
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nodes (i, j)). A cut, 8(S), of a set of nodes S is the set of edges with precisely one 
endpoint in S. A weight of &(S), cut(S), is defined to be ^^ =] Tu . The cut of a ball is 

the cut induced by the set of vertices included in the ball. The volume, vol(S), of a set of 
nodes S is the weighted distance of the edges with both endpoints in S, i.e. j^Tijdij . 

5 Finally, the volume of a ball is the volume of b(i, r) including the fractional weighted 
distance of edges leaving b(i, r). In other words, if (j, k) is a cut edge of ball b(i, r) with 
je b{i y r) , then (j, k) contributes r Jk x (r - dy) weight to the volume of ball b(i> r). For 
reasons known to one skilled in the art, an initial volume (seed) I is also included to the 
volume of every ball (i.e. ball b(i, 0) has volume I). 

10 [0047] In one embodiment of the present invention, rounding is carried out by first 
generating a graph G(V, E), \v\ = n , having fractional assignment to the variables dy 
obtained from the linear program. Next, suppose the volume of the entire graph is 

F == *~^ f j€V dun ! , where the updating cost of the fractional solution is represented by 

C U F. If the initial volume of the balls is F/n, balls are iteratively grown around arbitrary 
15 nodes of the graph until the cost of the cut is at most cln(n+l) times the cost of the 
volume. A location area is then generated, consisting of the nodes in this ball. The 
nodes are then removed from the graph, and the process is repeated. 

[0048] Table 1 summarizes the steps used to round fractional values generated by the 
linear program into integer values. 
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Table 1. 


Round (G(V.E). fdq /): 

// Variable Initialization. 

// Main loop 
while 3 a node/e H 
S ^0 

//Grow ball 
repeat 

r <— r + A 

until cut weight (b(i, r)) < c\n(n + 1) • vo\(b(i, r)) 
L « — Lkj S 

end 


[0049] It should be noted that the order in which the nodes are considered is indeed 
arbitrary. That said, the present invention can also be applied to a heuristic that chooses 
an order In a particular manner, as follows. 

5 [0050] In a further embodiment of the present invention, suppose c is some constant 
which will be determined later, and S = min{(dy - r) : y £ b(i, r), (do - r) > 0} is the 
remaining distance to the nearest vertex (among those with distance greater than zero) 
outside the current ball, then the process will terminate in polynomial time with a 
solution L that satisfies a set of constraints, and have a cost which is not much more than 
10 the fractional volume F. 

[0051] It should be noted that the termination condition on region-growing 
guarantees an (9(log n) approximation to the updating component. 

[0052] Let a = c\n(n +1), then 

Update _Cost(L) = —Cu ^cut(b) 

2 balls b 
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-Cua^olib) 

2 balls b 


— Cucc 
2 


< -C«a(2F) 

< CaF 

5 where the second line follows from the fact that balls found by the algorithm are disjoint. 
Note also that C U F is precisely the updating cost of the fractional solution. 

[0053] It can also be seen that the balls have radius at most 1/c. This fact follows 
from the following known lemma. 

[0054] Lemma: For any vertex i and family of balls b(i, rj, the condition cut (b(i, 
10 r)) < c\n(n + 1) x vo\(b(i, r)) is achieved for some r < lie. 

[0055] Proof: We proceed by contradiction. We first set a = cln(rc + 1) . Next, the 
ball is grown continuously from r = 0 to r = lie, and suppose throughout this process, 
cut(b(i, r)) > axvol(b(i, r)) . The incremental change in the volume becomes: 

d(vol(b(i,r))) = tf(2>*<&+ £>(r -</</)) 

j,keb j<=b,keb 

15 X d(rjk(r-dij)) 

je.b,k<zb 

X rjkdr 

j^bMb 

cut(b(I, r))dr 
= avol(b(i, r))dr 
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by assumption. The initial volume of a ball is, by definition, F/n, and the final volume is 
at most F + F/n if the ball covers the entire graph. Therefore, 

r F + F/n 1 / * ( v X x Ale 

Inn voK&JTr)) ^ 

and so ln(n + 1) > —a = ln(n + 1) . 

c 

5 [0056] In an additional embodiment of the present invention, the rounded paging 
cost is bounded by: 

PagejCost(L) = AC P £ 

balls b ijeb 


= -h-*c> Z EM'*)"* 

c ^ balls b ijeb 

where AjCp^^ J* m . ^ t (l-2/c)wj is a lower bound on the paging cost of a fractional 

10 because the radius of the balls is at most lie and, therefore, by the triangle inequity l-dy 
> 1-2/c for any nodes i and j that belong to the same ball. In an additional embodiment 

c 

of the present invention, this implies that our solution applies a approximation to 

c - 2 

the paging cost. 

[0057] In yet a further embodiment of the present invention, a final approximation 
15 factor is the maximum approximation factor of the two components, namely: 

max(cln(n + 1 ), ) = <9(log n) (9) 

c — 2 

[0058] It can be shown that Equation (9) satisfies size constraints as well. More 
specifically, Equation (9) represents a pseudo-approximation. A pseudo-approximation 
gives an approximate solution to a problem with slightly different parameters. In this 
20 case, the size bound parameters, and W fnax , are perturbed slightly. Specifically, a set 
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c 

of location areas are generated such that each location area has size at most K niax 

c — 2 

c 

and weight at most W niax , 

c-2 

[0059] Proof of the first of these statements is as follows. (Proof of the second is 
similar.) Because we know that V/eV: X/ev ^ ~ ^ -^ max * Fix i » therefore, 


c 


c~2 c-2j^ v 


EC 1 "*) 


> 


C - 2 ; G LA(i) 


— s h-- 


Z 1 

[0060] The last term is the size of LA(7). Note that c is an arbitrary constant. The 
10 larger c is, the closer the solutions are to true size bounds. However, if the 
approximation factor grows like c\n(n +1), overall costs may increase. The tradeoff is 
between being close to the size constraint on the one hand and having a better overall 
cost on the other. It should also be noted that if the bounds K max and W nxax are not 
specified (i.e. can be arbitrarily large), then the techniques presented above can be 
15 viewed as exact approximations. 

[0061] In still further embodiments of the present invention, rounding can be carried 
out with more general size and weight constraints. For example, constraints for each 
vertex / e V on the size and weight of a location area containing i, i.e., |LA(i)|<& and 

^ijeiA{i) Wj ~ Wi may ^ e applied. This requires adding many constraints exponentially to 

20 a linear program for each node i which force subsets containing i that exceed Ki and W, 
bounds to have some large The rounding techniques set forth above guarantee that no 
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location area has a large diameter. Therefore, it can be assumed that large subsets will be 
subdivided. 

[0062] Though the above discussion has focused on a few rounding techniques, it 
should be understood that others may be used without departing from the spirit or scope 
5 of the present invention. 

[0063] The above discussion has also assumed that cells to be grouped, partitioned or 
clustered do not form a special cluster referred to as a "line" graph. 

[0064] If, however, cells can be formed into a line cluster, then the present invention 
is also capable of generating optimal or minimum costs associated with these special 
10 types of cell clusters. 

[0065] For example, many times wireless networks are constructed along interstate 
highways or the like. It turns out that many of the cells form a line cluster. 
Economically, it is sometimes better to treat such clusters separately. In an alternative 
embodiment of the present invention, the network costs associated with a line cluster can 
15 be given by: 


[0066] The generation of costs using Equation (10) is based on the following 
assumptions and explanation. 

[0067] A line is formally defined as a graph where exactly two nodes have "degree 
20 1", where the degree of a node is the number of edges that have an end-point at that node, 
and the remaining nodes have degree 2. Consider a line G(V y E), connectivity matrix B 
and size bounds and W niax . Assume further that line nodes are indexed adjacently in 
increasing order from 1 to n, i.e., for every i, j, > 0 and by ± 0 only if \i - j\ = 1. The 
cost of the optimal LAP on the line graph G, induced by the first i nodes, [1../] e V is 
25 denoted by Cost( i), where Cost(0) = 0 is the cost of an empty line. 


i 



(10) 
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[0068] Cost(i) (and the actual LAP itself) can be computed recursively. For example, 
suppose the optimal LAP in a graph G, is {S It ... f S k } for some k. Let q be the index of the 
left border node in S* (i.e., the node with the lowest index in S*)< Then the optimal LAP 
for G q .j must have the same cost as {S/,...,S*} or else {Sy,...,S*} would not be optimal for 
5 Gj. Thus, Cost(i) is the sum of Cost(q-l) and the cost incurred by S* is that given by 
Equation (9). 

[0069] In yet another embodiment of the present invention, cells in a line may be 
grouped using a dynamic program representing a sum of weighted values associated with 
each cell and each edge between adjacent cells, and grouping constraints. Cells are then 
10 assigned to a group based on solutions of the dynamic program. Dynamic programming 
techniques are known by those skilled in the art. 

[0070] The dynamic program provided by the present invention provides an optimal 
LAP for a line graph in time 0(n ), where n is the number of nodes in the line. First, the 
Cost(i) for i from 1 to n is calculated and the result of each iteration is stored in a table so 
15 it may be used in the next iteration. In addition, the sum of weights is also stored to 
avoid an extra factor of n in run time. 

[0071] As before, connectivity and size constraints must be taken into consideration. 
Connectivity constraints by = 1 split the line graph into multiple line graphs. 
Connectivity constraints by = -1 simply forces nodes to be in same location area. 

20 [0072] Size constraints are easily accommodated as well. As an example, consider 
the maximal LA size constraint K niax (the weight constraint W max is analogous). When 
a = max{l,z - K max + 1}, the size constraint forces the border node q to be amongst nodes 
a,...,/, and so the minimization in Equation (10) is just taken over this range. 

[0073] Table 2. summarizes the steps involved in generating location area groups for 
25 a line ignoring connectivity constraints. 
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Table 2. 


//Variable Initialization. 
C[0] = 0 
L[0] = 0 
foj = 0 

//Main loop from 1 to n. 
for / = 1 to n do 
C777 = ~ 

//Updating the max weight of LA LA(i). 
a - max { 1, i - ^„ 2 ^ + 7/ 

//Loop for checking all border nodes. 
for <jr = a to i do 

//Checking if q can be a border node. 
if b q .} y q = 0 then 

tmp = A'C P -(i-q + l)'WiA + 

+ C w -/,-i,, + C[9-l] 
if rm/? <C/"/7 then 

//A less expensive LAP was found. 
C[i] = tmp 

L[i\=L[i]=L[q-l]^f{q..iJJ 
end-if 
end if 

Wia = W IA - W q 

end-for 
end for 

return C[n], L[n] 
end 


where L and C are two arrays that store the optimal LAP and its cost for the segment 
[1../], respectively and where the variable Wla records the total weight of the nodes q,...,i 

initialized by Wia = ^7 w/ , and decreased by w q before increasing a border node index, 
q. 

[0074] Given a line G(V, E), connectivity matrix B and maximal LA size Kmax, an 
optimal LAP of G(V, E) and its cost can be determined for a line. 

[0075] For an empty line without nodes there are no paging or updating costs and its 
total cost is Cost[0] = 0. In an additional embodiment of the invention, the cost of 
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different solutions is calculated when each one of the feasible border nodes q of LA(i) 
(the last max { 1, /- K niax + 1 } nodes of the line), is checked using the equation, 

A>C P -(i-q + l)-WiA + Cu'r<j-i, q + C[q-l] (11) 

[0076] By inductive assumption, C[q] y q < r, can be shown to be the cost of the 
5 optimal LAP for the line graph induced by nodes [l..<j]. Therefore, the cost of the 
optimal LAP can be generated when a given node q is forced to be a border node of 
LA(i). Because the minimum of these values is used over all feasible border nodes q, 
both the optimal LAP and its cost can be generated. 

[0077] The above discussion has presented some examples of the present invention. 
10 Further examples may be envisioned by those skilled in the art. For example, one or 
more of the functions, steps or processes discussed above may be carried out by one or 
more programmed devices (e.g., computer, microprocessor, etc. and/or associated 
memory devices, hard drives, floppy discs, etc. using software, firmware, etc.). The 
scope of the present invention is, however, defined by the claims which follow. 


